September 29, 2011 0:8 



WSPC - Trim Size: 9.75in x 6.5in qgp 



1 



^ I Fluid QCD approach for quark-gluon plasma in stellar 

^ ■ structure 

jyy ■ T.P. Djun" and L.T. Handoko"''' 

Group for Theoretical and Computational Physics, Research Center for Physics, Indonesian 
Institute of Sciences, Kompleks Puspiptek Serpong, Tangerang 15310, Indonesia 



00 
(N 



X 



^'^ Department of Physics, University of Indonesia, Kampus UI Depok, Depok 16424, Indonesia 



I ' The quark-gluon plasma in stellar structure is investigated using the fluid-like QCD 

^ l' approach. The classical energy momentum tensor relevant for high energy and hot plasma 

1^ having the nature of fluid bulk of gluon sea is calculated within the model. The transition 

. of gluon field from point particle field inside stable hadrons to relativistic fluid field in 

hot plasma and vice versa is briefly discussed. The results are applied to construct the 
equation of state using the Tolman-Oppenheimer-Volkoff equation to describe the hot 
. plasma dominated stellar structure. 
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o 

■ 1. Introduction 

o 

Recent experiments in the last decades on relativistic nuclear collisions shed light 
on the phenomena of hot plasma formed by dense quarks and gluons. Those exper- 
iments suggest that the quark gluon matter behaves more like a deconfined quark- 
gluon plasma (QGP) liquid. ^'^ A comprehensive review on this matter is given by 
E. Shuryak.3 

I This fact immediately encourages some models based on the (relativistic) hydro- 

dynamic approaches. In particular, dissipative ideal hydrodynamics has been used 
to fit some experimental data at high energy heavy ion program at the Relativistic 
Heavy Ion Collider (RHIC).^ The successful fit requires the models to take into 
account very small value of the ratio shear viscosity over entropy. However the 
puzzle must still be confirmed by the next coming experiments at the Large Hadron 
Collider (LHC)." 

Since QGP is containing many quark-anti-quarks and gluons, it is considerable to 
treat it using the well-established quantum chromodynamics (QCD). In pure QCD, 
QGP is described as a quark soup before hadronization which is a phase of QCD, 
and exists at extremely high temperature and/or density. It is argued that this phase 
consists of almost free quarks and gluons. Therefore, the phase transition from the 
deconfined QGP to the hadronic matters or vice versa gets particular interest in this 
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approach. It unfortunately turns out to the many body problems with large color 
charge which cannot be calculated analytically using perturbation. As a result, the 
main theoretical tools to explore the QGP within QCD is lattice gauge theory. The 
lattice calculation predicts that the phase transition occurs at approximately 175 

MeV.12.13 

On the other hand concerning that the QGP is a strongly interacting elemen- 
tary particle system which should be governed by strong interaction, while it also 
dissolves into an almost perfect dense fluid of quarks and gluons,^* it is plausible 
to describe it as a fluid system. In this sense, there are approaches based on unify- 
ing or hybridizing the charge field with flow fleld.^^"^^ Recently, some works have 
constructed the models in a lagrangian with certain non-Abelian gauge symmetry 
to the matter inside the fluid. ^"^'^^ 

In this paper, the energy momentum tensor of QGP within the recent fluid QCD 
model^^ is investigated. Further, the equation of state relevant for hot plasma dom- 
inated stellar structure is constructed using the so-called Tolman-Oppenheimer- 
Volkoff (TOV) equation. 

The paper is organized as follows. First we briefly introduce the underlying 
model of gauge invariant fluid lagrangian and discuss the relevant physical scale 
and region within the model. Then, the energy momentum tensor in the model is 
derived and investigated. Subsequently it is followed with relevant equation of state 
in a particular geometry using TOV equation.^^'^^ Finally, the paper is ended with 
a summary. 

2. The model 

Let us adopt the model developed by Sulaiman et.al. It describes the QGP as a 

strongly interacting gluon sea with the matters of quarks and anti-quarks inside. The 
model deploys the conventional QCD lagrangian with SU(3) color gauge symmetry, 
that is, 

£ = tQrd.Q - mgQQ - ^S^.S-"" + gs^.U'^'' ■ (1) 

Here Q and 11^ represent the quark (color) triplet and gauge vector field, gs is the 
strong couphng constant, = QT°-^nQ and T^'s belong to the SU(3) Gell-Mann 
matrices. The strength tensor is 5;^^ = d^U^ - d^V^ + gsf^'^UlU^ with is the 
structure constant of SU(3) group respectively. It should be noted that the quarks 
and anti-quarks feel the electromagnetic force due to the U(l) field A^, but the size 
is suppressed by a factor oi e/g = ^Jajas ^ O(10~^). 

Following the original model, the gluon fluid is put to have a particular form 
in term of relativistic velocity as,'^ 

c/; = (^7o^u'^)^w^</., (2) 



''This form was first proposed in tlic early work of Sulaiman et.al. in 2005 available in 
arXiv;physics/0508219. The form was then adopted in the work of Bambah et.al. in arXiv:hep- 
th/0605204, but the citation to the original work disappeared in their published version.^^ 
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lij^ = 7v"(l, v") and ^'^a = (1 — |v'^p) . (/) is a dimension one scalar field to 
keep correct dimension and should represent the field distribution. It is argued that 
taking this form leads to the equation of motion (EOM) for a single gluon field as 
follow,24 

^ (7v<.v«,^) + V (7v^</') = -9s j dx {JS + FS) , (3) 

where J^^ is the covariant current of gluon field, and is an auxiliary function 
which can be found in the original paper. It has been concluded that Eq. (3) 
should be a general relativistic fluid equation, since at the non-relativistic limit Eq. 
(3) coincides to the classical Euler equation. 

More precisely, Eq. (3) provides a clue that a single gluonic field may behave 
as a fluid at certain scale, beside its point particle properties with a polarization 
vector in the form of = (j). One can consider that there is a kind of "phase 
transition" , 

hadronic state i — > QGP state . (4) 

As the gluon field behaves as a point particle, it is in a stable hadronic state and is 
characterized by its polarization vector. On the other hand in the pre-hadronic state 
(before hadronization) like hot QGP, the gluon fleld behaves as a highly energized 
flow particle and the properties are dominated by its relativistic velocity. 

One should also recall that the wave function for a free particle satisfies 
_l_ jji^'^ _ Q'^Qi^]^ = Q with a solution ~ exp(— ipj^a;'') where is 
the 4- momentum. For a massive vector particle, i.e. mu 7^ 0, we have no choice 
but to take d^U/j, = 0. It is not a gauge condition like the case of massless particle. 
This then demands p^e^ = 0. Therefore the number of independent polarization 
vectors is reduced from four to three in a covariant fashion. In contrast with this, 
in the case of massless bosonic particles like gluon there axe only two degrees of 
freedom remain. Therefore, one should keep in mind that in the present model the 
spatial velocity has only 2 degrees of freedom, that means one component must be 
described by another two vector components. Fortunately, in real applications in 
cosmology or compact star, this requirement is satisfied by the assumption that the 
system under consideration is isotropic. 

From now, throughout the paper let us focus only on the gluon sea of plasma. 
This means one should consider only the related gluonic terms in Eq. (1), 

c, = -\s^^^s^'''' + g,j;u^r (5) 

3. Energy momentum tensor 

Now we are ready to proceed with deriving the energy momentum tensor within 
the model. It should be pointed out that once the hot (high energy) QGP state is 
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achieved, the system is assumed to be predominated by the classical motion rather 
than the quantum effects. 

Therefore the total action of matter for non-gravitational fields in a general 

geometry of space-time TZ is Sg = d'*x y/—gCg, where g is the determinant of 
metric g^,^. It is well-known that the variation of Sg in the metric is given by 
SSg = — 5 /tj d'^x ^/^Tfj.v Sg'^" . Since the energy momentum tensor density is, 

^Tg^gi^"' ^' 

one obtains, 

"TiJ-v — ^lip^'^u ~ Qui^^g + J^J/";/ • (7) 

It is clear that Eq. (7) is symmetric as expected to fulfill the Einstein gravitational 
EOM. The total energy momentum tensor T^j, is given by integrating out Eq. (7) in 
term of total volume in the space-time under consideration. T^i, is a result of bulk 
of gluons flow in the system. 

Furthermore, in a general space-time coordinates, the components of energy 
momentum tensor determine the total energy density (Tqo), the heat conduction 
{Toi,io), the isotropic pressure {Tu) and the viscous stresses (Tjj with i ^ j) of the 
gluonic plasma. Of course, in this case the derivative 9^ inside the strength tensor 
Sij,i, should be replaced by the covariant one, V^. Also, the energy momentum tensor 
satisfies the conservation condition, T'^" = 0. Nevertheless, one can trivially 
conclude that the model induces non-zero viscosity since generally Tjj ^ for i ^ j. 
From the experimental clues, however the size should be small such that it is always 
treated perturbatively in most hydrodynamics models. ^"^"^ 

Before going further to apply these results, one should determine the quark 
current J° in Eq. (7). This can be simply calculated by considering the EOM (Dirac 
equation) of a single colored quark (q) or anti-quark (q) with 4-momentum p^. Since 
the solution of the EOM is q{p, x) = u{p) exp{—ip-x), one immediately gets U7^m = 
4p^. Assuming that all colored quarks / anti-quarks have the same momenta and 
the velocity of gluons are homogeneity, approximately J^J/"'* oc 4p^{7'* = Amqcj) 
since m^m'' = = 1. 



4. Equation of state for stellcir structure 

This section is dedicated to provide an example on the applications of the present 
model to describe the stellar interior, in particular the compact stars which are still 
dominated by hot plasma before transforming itself into neutron star. 

The stellar structure is commonly described as a static spherically symmetric 
space-time represented by Schwarzschild geometry. This means one deals with the 
relativistic gravitational equations for the interior of spherically symmetric plasma 
distribution. In the region under consideration the presence of the flow gluonic flelds 
induces non-zero energy momentum tensor which is making up the star. This is the 
phase before the neutron star is getting mature. Starting from the stellar nebula 
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made of hot plasma which is gradually getting colder as the hadronization occurs 
from the colder surface, while the inner core is still in pure hot QGP state. 

As a consequence of the diagonal metric of Schwarzschild space-time, the model 

falls back to the perfect fluid without viscosity and heat conduction, i.e. T^i = 
Tij = for i ^ j. Also, since the plasma distribution should be spherically isotropic, 
it is considerable to put vi = V'2 = V3 = v as constant for all colored gluons. 
This assumption is consistent with the degree of freedom counting discussed in 
the preceding section. Moreover, the vanishing off-diagonal components of the Ricci 
tensor, i?io, actually forces the spatial 3-velocity of the fluid must vanish everywhere. 
Hence particular assumption for Vi is indeed not necessary. However, the gluon 
distribution still depends on the radius length, (j) = 4>{r). 

For the sake of simplicity one can put homogeneous gluon fields for all color 
states, i.e. Uf^ =U^ for all a = 1, • • • ,8. This yields. 



1 



^9sfQmQ(l>{r)--alfl(t>{rf 



9nv ) 



(8) 



where fg is the factor of summed colored gluon states from the structure constant 
jobc^ while /q is the factor of summed colored quark states from J'^U"'^ . Remind 
that the energy momentum tensor for perfect fluid takes the form. 



Tfii, = {£ + V) u^u^ -V g^v 



(9) 



Here £ and V denote the density and isotropic pressure for single fluid field, each 
is related to the total density and pressure of the system through p = § d^x £ and 
P = § di^xV respectively. Obviously, from Eqs. (8) and (9) one can obtain the 
density and pressure in the model as follows, 



P{r) 



dt / dV 



4:gsfQmQ(l){r) 



1 



^gsfqmQ 



dV 



9 s fg 

16 /q mq ' 



(r), 



(10) 



p{r) 



dV 



^9sfQmQ<l>{r) + -gif^cl>{rr 



■^9sfQmQ 



I 



dV 



^9sfg 

16 /q mg 



(11) 



at a finite temperature (3 = 1/T in a 3-dimensional spatial volume V. 

The proper spatial volume element for Schwarzschild geometry is dV = 
\/ B{r)r^ sin 9 dr d9 dip with radius r and two angles 9 and if in spherical coor- 
dinates. The solution for B{r) is given by. 



B{r) = 



1 - 



2GTO(r)" 



(12) 



September 29, 2011 0:8 



WSPC - Trim Size: 9.75in x 6. Sin qgp 



and m(r) = 47r dfp{f) is the 'bare mass'. This generates the proper integrated 
mass m(r) contained within a coordinate radius r inside the star. On the other 
hand, the stellar structure with Schwarzschild geometry is well known as the TOV 
equation which relates density and pressure in a unique way,^^'^^ 



dP(r) 1 



P(r)] [47rGP(r)r^ +Gm(r)] 



1 



2Gm{r) 



(13) 



Substituting Eqs. (11) and (10) into Eq. (13) provides a direct relationship be- 
tween density and pressure. In another word, one can obtain a contour of the equa- 
tion of state of gluonic plasma in term of distribution function <^{r) and temperature 
T. 



5. Summary 

The phase transition between stable hadronic state and highly energized QGP state 
is briefly discussed using recently developed fluid QCD lagrangian. It is argued that 
inside the hadronic state the gluons behave as point particles and the properties are 
determined by its polarization vectors. However, in the hot QGP state, the gluon 
should behave as fluid particles and characterized by its relativistic velocities. 

Then, the energy momentum tensor for the fluid particle is investigated. In 
particular, a detailed derivation has been worked out in the case of Schwarzschild 
space-time that is relevant for stellar structure of unmatured stars composed of hot 
QGP. In principle one can obtain a kind of equation of state of gluonic plasma inside 
the star through the TOV equation. This approach would enable us to describe 
the stellar structure without assuming a particular relation between pressure and 
density. 

Further works, especially on the numerical analysis, are still in progress and will 
be reported elsewhere. 
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